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Correlation Effect on Peierls Transition 
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The effect of correlation on Peierls transition, which is accompanied by a dimerization, td, of a 
bond alternation for transfer energy, has been examined for a half-filled one-dimensional electron 
system with on-site repulsive interaction ([/). By applying the renormalization group method to 
the interaction of the bosonized Hamiltonian, the dimerization has been calculated variationally 
and self-consistently with a fixed electron-phonon coupling constant (A) and it is shown that td 
takes a maximum as a function of U . The result is examined in terms of charge gap and spin gap 
and is compared with that of the numerical simulation by Hirsch [Phys. Rev. Lett 51 (1983) 
296]. Relevance to the spin Peierls transition in organic conductors is discussed. 
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Peierls transition for a one-dimensional half-filled elec- 
tron system coupled with phonon has been studied exten- 
sively since the Su-Schricffcr-Hccger model was proposed 
for the quasi-one-dimensional conductor polyacetylene 
with a bond alternation^' The role of correlation in 
such a model has been examined by introducing an on- 
site repulsive interaction, U . The calculation using the 
Hartree-Fock approximation leads to the Peierls state, 
which exists only for U smaller than a critical value of 
the order of the band width.B'B This study has been fur- 
ther developed by taking into account a one-dimensional 
quantum fluctuation. The numerical simulation exhibits 
an enhancement oLdimta-i^jation in the presence of re- 
pulsive interactionD'BQ'lJ'ElLF A notable finding is that 
the dimerization takes a majdmum at a value of U be- 
ing nearly the band widthErBtl'' and that the charge gap 
becomes much larger than the dimerization gap.© Re- 
garding the case of weak coupling, the analytical method 
of a renormalization group (RG) has also exhibited the 
enhancement of td as a function of Ux3 The effect of 
finite phonon frequency has been exploLCiLusing the RG 
method based on the bosonization.OEj'Ej-' It has been 
shown that the half-filled case leads to a competition 
between the state with both spin and charge gaps and 
the state with only a charge gap on the plane of the 
phonon. frequency and the electron-phonon coupling con- 
stant .113 In contrast, the case of strong coupling with 
large U has been examined in terms of the opposite ap- 
proach, i.e., the expansion of \/U) and mapping the elec- 
tron system into a spin 1/2 chain system, which leads to 
the spin Peierls transition. The successful treatment of 
the quantum fiuctuation leads to the occurrence of the 
spin Peierls transition for an arbitraxij-Lnagnitude of the 
electron-phonon coupling constant.llJ'lla'y However, it 
is not clear why the optimum condition for the Peierls 
transition is given by the intermediate coupling of U . 

In the present paper, the unconventional role of U on 
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the Peierls state is examined by applying the RG method 
to the bosonized Hamiltonian. It is demonstrated that a 
maximum of the dimerization occurs for /7 ~ 2t {At is the 
band width) for weak electron-phonon coupling (A), and 
that the charge gap is well separated from the spin gap at 
the maximum. Furthermore, we discuss the relevance of 
this to the spin Peierlaixansition in an organic conductor, 
(TMTTF)2PF6 sah.By) 
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where ^ denotes a creation operator of a conduction 



electron with spin a{=1, |) at the lattice site j. The first 
term is the kinetic energy, where t is the uniform transfer 
energy. The quantity, td, which denotes the dimerization 
due to Peierls distortion, is determined so as to minimize 
the total energy. The U term with rij^a- = cj ^Cj,a denotes 
an on-site repulsive interaction, and the last term with 
a constant C is the elastic energy for_Lha.distortion. 

Applying the bosonization methodJlj'Ej' eq.(l) for the 
half-filled case is rewritten as 
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yw sin 6*+ cos(/)+ + (l/8A)j/^ 
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where 0±{x) and (/>±(a;) denote phase variablfis for 
the charge and spin fiuctuation, respectivelyJM' and 
[9+{x),d^{x')] = [0+(x),0_(x')] = i7rsgn(x-x') . 
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The 2/1/2 term and the yw term represent the Umk- 
lapp scattering and the dimerization, respectively. The 
coefficients in eq.(2) are given by Kp = (1 + C/)~^/^ 



\l/2 



= {l-U)-y\ Vp = vf{1 + C/)i/2, = fF(l - U) , 
yi/2 ^ Va ^ U, yw ^ 8atd/vF, 1/A = TrvpC/Sa, vf = 
2tasiiikFa(— 2ta), /cp = 7r/{2a) and U — Ua/ (ttvf)- 
The quantity a is a cutoff parameter of the order of the 
lattice constant, a; the ratio of a to a will be evaluated 
later. The quantity y^ is determined by minimizing the 
total energy, and the resultant self-consistency equation 
is written as 

^ yvy = ( sin 6*+ cos 0+ ) (= A) , (3) 

where yw — 'itda/{ta). The quantity A, which denotes 
an order parameter for the dimerization, is calculated 
self-consistently in terms of eq.(2). Equation (2) is exam- 
ined using the RG method with a scaling a — > a{l + dl). 
Within the lowest order of perturbation^ pRQ^cqjiMkms 
for the coupling constants are derived asE§E3llil'EJ'E3c^ 
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t; yi/2{l)yw{l) - \ya{l)yw{l), 



(4) 
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where Ga{l) — 2{Ka{l) — 1) and initial values at I = 
are given by coefficients of eq.(2). The relevance of yi/2, 
2/o- and yw corresponds to the appearance of charge gap, 
spin gap and dimerization, respectively. 

For calculating the r.h.s. of eq.(3), we use the fact that 
the derivative of (sin cos 0+) with respect to yw can 
be expressed in terms of the response function given by 



R{\ri ^ r2\) ^ {TrOd{ri)Od{r2)) , 



(5) 
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Fig. 1. Z-dependence of R{1) of eq.(6) for = 0.1, and U =0, 
1, 2, 3 and 4 where R(l) decreases and takes a minimum at 
I = Im- t = 1 and a/a = 1.52. The arrow denotes I = 1^ 
given by y\/2{lc) = 1- In the inset, the spin Peierls (SP) region 
and the Peierls (P) region correspond to 1^ < Im and Ic > im, 
respectively, where the boundary denotes a crossover between 
these two regions. 



where I = ln([a;^ 4- [vpTy]^/"^ /a). In eq.(6), we intro- 
duced I = Ic with 2/1/2(^0) = 1, from which the charge 
gap is calculated. The response function for I > Ic is 
calculated by estimating (sin^^.) in the presence of the 
charge gap. Although one expects i?(oo) ^ A^, we use 
the following method to estimate A within the present 
RG equation, which leads to R{1) taking a minimum at 
I = Im- Since R{1) as a function of I becomes invalid for 
/ > Im, we use a self-consistency equation for A(= A{td)) 
given by 



Mtd) = / dtl 



td 



dl'e'' {R{t',J') - A2(0}(7) 



whereOd(r)(= sin 6*+ (f ) cos (r) — — (7ra/4a) J2cri~^y^] 
h.c.) denotes an operator for dimerization. The vector 
r consists of space x{— ja) and imaginary time r while 
T-T is the time-ordering operator. _LIsing the RG method, 
R{td,l){= R{\r\)) is calculated as@ 



R{td, = 7 exp 



dl' {Kp{l') 



+K„{l')~yi/2{l') + ya{l')) , {l<lc) 



R{td,l) 



- exp 



dl' {Kp{l') 



+K„{l')^yi,2{l') + y.{l')) 



X exp 



dl'{Ka{l')+yM , {l>lc) 

(6) 



We assume t = 1 in the following numerical calcula- 
tion, if there is no confusion. For estimating the r.h.s. 
of eq.(7), we calculate RG eq.(4). When I increases with 
jix^djiJi{o^'0)i~one finds the decrease of Kp{l), G^il) and 
2/(t(0 and the increase of 2/1/2 and yw- Behaviors at large 
I are classified into two groups depending on the magni- 
tude of yw{Q)- For small 2/w(0), 2/1/2(0 is larger than 
yw{l)-, while 2/1/2(0 becomes smaller than yw{l) with 
increasing yw{^)- The former case denotes the Peierls 
state in the presence of a well-developed charge gap since 
the relevance of 2/1/2 {yw) corresponds to the formation 
of the charge gap (the formation of both charge and spin 
gaps due to the Peierls distortion). 

The response function is calculated by substituting the 
solution of RG equations into eq.(6). In Fig.l, the l- 
dependence of R{1) is shown for td = 0.1 and certain 
values of U . The results are shown in the region with 
I < Im, where R{lm) is a minimum of R{1). The long 
range order of the Peierls state is related to the rele- 
vance of yw, which leads to the appearance of Im in the 
present approximation of the RG equations. The arrow 
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for [/ = 2, 3 and 4 denotes the location at I = Ic where 
y 1/2(^0) — 1- From a comparison of Ic with Im, one finds 
a distinction between case I { Im > Ic ) and case II ( 
Im < Ic)- The charge gap is aheady developed prior to 
obtaining the Peierls state in case I. The inset depicts 
the spin Peierls region (case I) and the Peierls region 
(case II) on the plane of U and t^. The charge gap dom- 
inates for and td < 0.5. The solid curve shows 
the boundary as a crossover between these two regions, 
where the state moves continuously. 

An estimation of the ratio of a/ a is needed for the 
calculation of A of eq.(7) as a function of td and 
U. When U = 0, A becomes equal to the order 
parameter of the conventional Peierls distortion and 
thus the r.h.s. of eq.(3) can be calculated rigorously 
as A — > Apcioris = (tdOf/a) /J" dz (sin z)^/ [(2 COS z)2 + 
(2idsinz)^]^/^ , which leads to (tda/a) x ln(1.4/id) for 
small td- By comparing Apcieris with that obtained from 
eq.(7), we obtain a/a ~ 1.5 ± 0.02 for small td{ ^ 0.2). 

In terms of the response function, A is calculated self- 

0.2| ■ , ■ , ■ , ■ , ■ 1 
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Fig. 2. 4d-dependence of A calculated self-consistently by using 
eq.(7) for t/ = 0, 1, 2, 3 and 4. The dashed curve denotes 

^Peierls- 

consistently from eq. (7). In Fig. 2, A is shown as a 
function of td for certain values of U . The dashed curve 
denotes Apoicris which coincides well with that of J7 = 
(solid curve) for ^ 0.2. With increasing td, A for fixed 
U increases monotonically. The self-consistent result of 
A obtained from eq. (7) is smaller than i?(/„i)^/^, e.g., 
A/i?(/™)i/2 = 0.43 ,0.45 ,0.52, 0.47 and 0.39 for U = 
,1, 2, 3 and 4. This is reasonable since the lowest order 
RG equation with the relevant coupling overestimates 
the response function for large I. From the comparison 
of A for [7 = 2 with that for U = 3, it is found that A as 
a function of t/ ( ^ 2) decreases, indicating the existence 
of an optimum value of U for the Peierls state. 

Based on Fig. 2, we obtain the dimerization, td, repre- 
senting the Peierls state, which is calculated using eq.(3). 
In Fig. 3, is shown as a function of A with certain val- 
ues of U. The dashed curve denotes td{— 1.4exp[— 1/A]) 




Fig. 3. A-dependence of td, which is calculated from eq.(3), for 
(7 = 0, 1, 2, 3 and 4 where the dashed curve denotes Apoicris 
corresponding to the conventional calculation in the absence of 
U. 



obtained from Apoicris: which coincides with td in the 
case of [/ = (solid curve) for 0.25^A^0.5 within the 
visible scale. For small A, td decreases exponentially but 
is enhanced sufficiently for U — 1, 2 and 3. We note 
that td takes a maximum around ?7 ~ 2. Such a U- 
dependence of td, indicating a crossover from the regime 
of weak coupling to that of strong coupling, is elucidated 
by examining the relevant term of yw sin 9^ cos 0+ in 
eq.(2) and the corresponding RG equation (the last equa- 
tion of eq.(4)). For small U, the yvv'-term is enhanced by 
the Umklapp scattering of the j/i/2-term, which leads to 
the suppression of charge fluctuation (9+). For large U, 
both j/cr and increase the quantum fluctuation for the 
spin part {(f)^ ) resulting in the suppression of < cos > 
for the Peierls state. 

The {/-dependence of td{= AXa/a) is obtained by ex- 
amining td with fixed A in Fig. 3. In Fig. 4, td with 
A — 0.25 is depicted as a function of U by the solid 
curve. As expected from Fig. 2 and Fig. 3, td takes a 
maximum around U = 2. The triangle, which represents 
the numerical result obtained by IIirsch,tP is compared 
with the dashed curve of the present calculation with 
A = 0.42 (dashed curve) corresponding to the parame- 
ter of the numerical simulation. The fact that the dashed 
curve is larger than the numerical result (triangle) is rea- 
sonable since the present calculation treats the dimeriza- 
tion ( the bond alternation) classically while the numeri- 
cal simulation has been performed in the presence of the 
quantum fluctuation. We note that there is not much 
variation of the location of U for the maximum, e.g., the 
maximum appears at U — 2.2±0.1 in the interval range 
of 0.2 ^A^ 0.5. In the inset, the JJ-dependences of the 
charge gap (Ap) and spin gap (A^.) for A — 0.25 are 
shown by a solid curve and dotted curve, respectively. 
These gaps are calculated by Ap/cuc = exp[— ^c] and 
Ac/lUc = exp[-l^] with \yaila)\ = 1, where a;c(— 2.37) 
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td 0.2 




Fig. 4. Quantity td as a function of U for A = 0.25 (solidi curve) . 
The symbol (triangle) is the result calculated by HirschfP which 
is compared with the present calculation with A = 0.42 . The 
inset shows the charge gap (Ap) and the spin gap (Ao-) for A = 
0.25. 



is chosen so as to obtain Ap with td 



to the weU-known exact one. With increasing U, 



nearly equal 
Ap in- 
creases monotonically while A^ takes a maximum. The 
maximum of td is found to be accompanied by the sepa- 
ration of charge gap and spin gap. Thus, it is expected 
that with increasing U, there is a crossover from the con- 
ventional Peierls state with Ap ~ A^- to the spin Peierls 
state with the charge gap being much larger than the 
spin gap. 

When U is extremely large, the effect of td can be 
examined-using the bosonized phase Hamiltonian of spin 
1/2 chainta' with the antiferromagnetic exchange energy, 
4:{t±td)'^ /U . The calculation oitd, similar to eq.(3), with 
only a spin degree of freedom shows that td decreases 
monotonically with increasing U. From comparison of 
the present result with this limiting one (not shown ex- 
plicitly), the present calculation for [/ ^ 4 in Fig. 3 seems 
to be reasonably extrapolated to that of the spin 1/2 
chain except for small td (and then small A), which re- 
quires much accuracy of numerical evaluation. 

Finally, we discuss the spin Peierls transition in an 
organic conductor, (TMTTF)2PF6, which has..hccn ob- 
served by magnetic and X-ray experiments. E3'E3^ This 
conductor, which indicales the charge gap being much 
larger than the spin gapEJ' is often analyzed in terms of 
a model with spin_L/2 chain although U is of the order 
of the band width.ES' The conductor has a quarter-filled 
band with a dimerization, which may be considered to 
be half-filling. However, such an effectively half-filled 
band strongly reduces the magnitude of Umklapp scat- 
tering, resulting in the suppression of the Peierls state. 
In fact, the Umklapp scattering is estimated to be yi^2 = 
2{xd/ {l+x^))U , where vp = V2to,a^d Xd corresponds to 



a dimerization for a quarte£.-filling.E3' For Xd — 0.1 corre- 
sponding to TMTTF salt,llIP we find that td as a function 



of U decreases for small U while separation between Ap 
and Act appears for [/ ^ 3. Moreover, with decreasing Xd 
and fixed A( ^ 0.5), td decreases strongly for U = 4 '--^ 6. 
When U = 5.6 and A — 0.25(0.5), we obtain the dimer- 
ization as td = 0.036 (0.23), 0.019 (0.17) and ~ (0.05); 
the charge gap as Ap = 0.31 (1.05), 0.12 (0.73) and ~ 
(0.17); the spin gap as A,, = 0.18 (0.77), 0.08 (0.57) 
and ~ (0.15), for Xd = 0.2, 0.1 and 0, respectively. 
Based on such a consideration, the spin Peierls state in 
the organic conductor could be realized when Umklapp 
scattering induced by the electronic correlation becomes 
large. 

In summary, we have examined the effect of on-site 
repulsive interaction, C7, on the Peierls state with dimer- 
ization td- For small (large) U, td increases (decreases) 
where the charge gap is almost equal to (much larger 
than) the spin gap. The maximum of td indicates a 
crossover from the weak coupling regime, in which the 
Umklapp scattering suppresses the charge fluctuation, 
into a strong coupling regime, in which the quantum spin 
fluctuation reduces the effect of dimerization. 
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